1 è®åºæ¬æªæ¡
# read files
pacman::p_load(readr, dplyr, sf, tmap, stringr,tidyr,ggplot2,plotly)
load("hen.rdata")
load("GeoTW108.rdata")
load("RiTW108.rdata")å ±äžåè³ææªïŒ
(1)hen.rdataïŒ2018é«éåžé·éžè & 2020瞜統倧éžïŒé«éåžåéä¹æç¥šçãéåžé·åŸç¥šçã
(2)GeoTW108.rdataïŒå瞣åžåååºæ¬äººå£çµ±èšè³æ(sf)ã
(3)Ritw108.rdataïŒåºæ¬äººå£çµ±èšè³æïŒå å«è®æžè¡šã
2 åéæžéçµ±èš
## .
## äžæ°å 倧瀟å 倧寮å 倧暹å å°æž¯å 仿Šå å
§éå å
éŸå
## 86 9 25 18 39 16 18 12
## å·Šçå æ°žå®å ç°å¯®å ç²ä»å ææå é£çªå€å å²¡å±±å æåå
## 39 6 10 7 8 3 34 24
## é¿è®å åéå åé®å çŸæ¿å èé
å èæå èè£å æ¡æºå
## 12 20 60 19 70 3 15 8
## æ¢å®å é³¥æŸå æ¹å
§å æ°èå æ¥ æ¢å 路竹å éŒå±±å æå±±å
## 15 7 14 32 37 20 39 21
## ææŽ¥å é³³å±±å æ©é å çå·¢å åœéå 鹜åå
## 14 76 17 11 12 22
3 æç¥šç & åŸç¥šççæ³(以åå±±ååäŸç)
M = c("鳳山å","æå±±å","岡山å","éŒå±±å")
K1 <- K %>% filter(TOWN %in% M) %>% .[,-1]
facet = c("Rate18","Rate20","Hen18","Hen20")
tm_shape(K1) + tm_polygons(facet,popup.vars = c("TOWN","VILLAGE"))+tm_facets(ncol = 2,nrow = 2)+
tmap_mode("view")## tmap mode set to interactive viewing
4 åå±±åééåžé·åŸç¥šå·®ç°(2018 VS 2020)
幟ä¹çºè² å·®ç°ïŒæ
å¯ä»¥ççåªäºå°åæ¹è®çæ³å€§æå°
tm_shape(K1)+
tm_polygons(col = "Delta",
popup.vars = c("VILLAGE","Delta"),palette = "Greens",style = "equal")+
tmap_mode("view")## tmap mode set to interactive viewing
5 æ¢èšæ¶èæ¯è18,20幎éåžé·åŸç¥šå·®ç°å
è²å¡çºéåžé·åŸç¥šçå·®ç°
ååæ·±æ·ºçºæ¶èæ¯ïŒè¶æ·±ä»£è¡šè¶é«
# æ¢èšæ¶èæ¯å°åè18,20幎éåžé·åŸç¥šå·®ç°å
tm_shape(K)+
tm_polygons(col = "Delta",
popup.vars = c("TOWN","VILLAGE","Delta"),palette = "Greens",style = "fixed",breaks = seq(-0.35,0,0.05))+
tm_shape(K)+
tm_symbols(col = "FLD3",border.col = "white",size = 0.005)+tmap_mode("view")K2 <- sf::st_transform(K, "+init=epsg:4326") %>% st_centroid() %>%
as.data.frame()%>%separate(geometry,c("lng","lat"),",") %>%
mutate(lng = as.numeric(gsub("[A-Za-z\\(]+","",lng,)),
lat = as.numeric(gsub("[\\)]","",lat))) %>% na.omit()## Warning in st_centroid.sf(.): st_centroid assumes attributes are constant
## over geometries of x
## Warning in st_centroid.sfc(st_geometry(x), of_largest_polygon =
## of_largest_polygon): st_centroid does not give correct centroids for
## longitude/latitude data
st_transform()ïŒèœæç¶ç·¯åºŠã
st_centroid()ïŒèšç®åæ¬polygonsçäžé»ã
6 æç¥šäººå£å100åçæç¥šè®åçæ³
æŸåºå¯æç¥šäººå£å100éé
䞊äžäŸç
§å898åæéå
©å¹ŽåŸç¥šå·®ç°çäžäœæžååçºå
©åé¡è²ïŒç¶ (è®å倧)ïŒè(è®åå°)
# æç¥šäººå£å100åçæç¥šè®åçæ³
pacman::p_load(spatstat,ggplot2,leaflet,dplyr,htmlwidgets,stringr,ggmap,leaflet.extras,viridisLite,htmltools,sf)
lat = median(K2$lat);lng = median(K2$lng)
K3<-K2 %>% arrange(-Pop20) %>% head(100)
picon = pulseIcons(color = ifelse(K3$Delta > quantile(K2$Delta,0.5),"blue","green"))
K3 %>% leaflet() %>%
addTiles() %>% addPulseMarkers(icon = picon,label= ~ paste0(TOWN,VILLAGE,paste0(round(Delta,4)*100,"%")))%>%
setView(lng,lat,10)7 æ¶èæ¯å100åçæç¥šè®åçæ³
# æ¶èæ¯å100åçæç¥šè®åçæ³
K3<-K2 %>% arrange(-FLD3) %>% head(100)
picon = pulseIcons(color = ifelse(K3$Delta < quantile(K2$Delta,0.5),"green","blue"))
K3 %>% leaflet() %>%
addTiles() %>% addPulseMarkers(icon = picon,label= ~ paste0(TOWN,VILLAGE,paste0(round(Delta,4)*100,"%")))%>%
setView(lng,lat,10)